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ABSTRACT 


/ 
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- — ^An  International  oceanographic  experiment,  Oonde  Va?,  Mas  conducted  In  1982  In 
the  Alboran  Sea  In  which  a  specialized  portion  was  dedicated  to  examining  the  spa¬ 
tial  and  temporal  variability  of  bio-optical  properties  around  the  Alboran  Sea  tyre. 
The  circulation  of  the  Alboran  Sea,  as  characterized  by  Atlantic  Inflow  (through  the 
Strait  of  Gibraltar),  the  Alboran  tyre,  and  coastal  water  masses,  was  analyzed 
through  use  of  Nimbus-7  Coastal  Zone  Color  Scanner  (CZCS)  Imagery,  NOAA-7  Advanced 
Very  High  Resolution  Radiometer  (AVHRR)  Imagery,  and  aircraft  and  ship  data.  Diffuse 
attenuation  coefficients  calculated  from  CZCS  data  are  In  agreement  with  the  ship 
data.  A  strong  correlation  between  the  surface  temperature  and  ocean  color  fronts 
was  observed  In  the  two  satellite  Imagery  sets. 

A  6-day  sequence  of  CZCS  Imagery  shows  that  significant  changes  In  the  diffuse 
attenuation  coefficient  and  phytoplankton  pigment  concentration  occurred  In  the 
frontal  regions.  These  changes  are  attributed  to  the  bio-optical  horizontal  and 
vertical  variations  that  occurred  within  the  first  attenuation  length. 
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BIO-OPTICAL  VARIABILITY  IN  THE  ALBORAN  SEA 
AS  ASSESSED  BY  NIMBUS-7  COASTAL  ZONE  COLOR  SCANNER 


I.  INTRODUCTION 

Donde  Va?,  an  International  oceanographic  experiment  conducted  In  the  period 
June  through  October  1982  In  the  Alboran  Sea,  had  as  Its  major  focus,  understanding 
the  regional  circulation  and  the  physical  processes  that  force  the  dynamic  changes 
(La  Vlolette  et  al.,  1982  and  Kinder,  1983).  Analysis  of  Nimbus-7  Coastal  Zone  Color 
Scanner  (CZCS)  data  form  a  specialized  portion  of  the  experiment,  which  Is  aimed  at 
determining  the  spatial  and  temporal  variability  of  bio-optical  properties 
associated  with  the  circulation  and  the  frontal  dynamics.  Preliminary  results  of 
this  effort  are  described  In  this  paper.  An  important  aspect  of  the  effort  Is  an 
assessment  of  the  utility  of  existing  CZCS  algorithms  to  obtain  quantitative 
measurements  of  the  optical  diffuse  attenuation  coefficient  and  the  phytoplankton 
concentration. 

The  West  Alboran  Basin  Is  the  first  Mediterranean  basin  east  of  the  Strait  of 
Gibraltar  (Figure  1).  A  300-m  sill  located  within  the  Strait  forms  the  western 
boundary  of  the  basin,  and  a  series  of  subsurface  ridges  form  the  eastern  boundary. 
Along  the  European  and  African  coastal  boundaries  little  shelf  area  exists,  with  the 
bathymetry  dropping  to  depths  of  1000  m  approximately  5  km  offshore. 

The  general  circulation  In  the  Alboran  Sea  has  been  characterized  by  several 
Investigators  (e.g.,  Lannolx,  1974;  Parllla  and  Kinder,  1984;  and  La  Vlolette, 

1984).  They  describe  the  Atlantic  water  flowing  Into  the  Alboran  Basin  through  the 
Strait  of  Gibraltar  as  essentially  being  confined  to  the  upper  200  m.  This  water 
flows  eastward  along  the  Spanish  coast  for  approximately  100  km,  then  turns  south¬ 
ward  toward  a  point  slightly  west  of  Cap  Tres  Forcas  on  the  north  African  coast 
where  It  splits  Into  a  east-west  flow.  This  circulation  creates  a  large  antlcyclonlc 
gyre  (aptly  called  the  Alboran  Sea  tyre)  that  often  occupies  most  of  the  West 
Alboran  Basin. 

The  coastal  water  along  the  southern  Spanish  coast  Is  cold,  high  In  salinity, 
and  rich  In  nutrients.  Its  complex  nature  Is  Influenced  by  the  movements  of  the 
Atlantic  Inflow  water  as  well  as  the  local  meteorology,  bathymetry,  and  tides.  La 
Vlolette  (1984),  using  sea  surface  temperature  values  and  patterns  from  thermal 
Infrared  satellite  data,  showed  the  Atlantic  Inflow  water  has  associated  mesoscale 
cold-water  thermal  features  that  appear  to  be  cyclic  responses  to  the  tides.  La 
Vlolette  and  Kerllng  (1983),  using  aircraft  expendable  bathythermographs  (XBTs), 
showed  that  the  geographic  location  of  the  Atlantic  Inflow  Is  highly  variable.  The 
thermal  characterization  of  the  front  formed  between  the  waters  of  the  cold  Atlantic 
Inflow  and  the  warm  Alboran  Sea  tyre  Is  easily  distinguishable  by  the  gray  shading 
of  the  thermal  Infrared  satellite  Imagery.  However,  the  demarcation  between  the 
northern  edge  of  the  Atlantic  Inflow  and  the  coastal  water  Is  difficult  to  determine 
since  similar  cold  surface  temperatures  are  present  In  both  areas.  The  discrimina¬ 
tion  of  the  Coastal,  Atlantic  Inflow,  and  Alboran  Sea  tyre  waters  should  be  Improved 
by  combining  the  analysis  of  color  and  thermal  signatures  of  the  water  masses. 

The  larae  sea  surface  temperature  variability  such  as  described  by  La  Vlolette 
and  Kerllng  (1983)  and  La  Vlolette  (1984)  Indicates  strong  mixing  or  upwelllng  that 
should  produce  rapid  changes  In  the  amount  of  nutrients  available  for  biological 
activity.  In  addition,  the  volume  of  the  resulting  biomass  would  depend  largely  on 
the  time  available  for  the  organisms  to  respond  to  the  advectlon  or  depletion  of 
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nutrients.  Since  water  optical  properties  relate  directly  to  the  concentration  of 
phytoplankton,  their  bloom  (or  decay)  significantly  changes  the  optical  character¬ 
istics  of  the  water.  Thus,  sequential  CZCS  lamgery  should  reveal  changes  In  the 
surface  bio-optical  properties  In  a  fashion  similar  to  the  temperature  changes 
observable  In  sequential  thermal  IR  Imagery. 

Before  the  data  and  their  analyses  are  presented,  certain  basic  premises  must 
be  emphasized.  The  upwelled  radiance  (or  color)  of  the  ocean  results  from  back- 
scatter  within  the  water  column  combined  with  spectral  absorption.  Most  Important 
for  these  discussions,  approximately  90  percent  of  the  spectral  water-leaving  radi¬ 
ance  at  the  surface  Is  backscatter  from  within  the  first  optical  attenuation  length 
(Gordon  and  McCluny,  1975).  This  length  Is  defined  to  be  the  penetration  depth  at 
which  the  radiance  falls  to  1/e  of  Its  value  just  beneath  the  surface.  For  practical 
purposes,  therefore,  the  upwelllng  spectral  radiance  at  the  water  surface  sensed  by 
CZCS  may  be  considered  to  be  a  result  of  the  Integrated  effects  of  backscatter  and 
absorption  within  the  upper  attenuation  length,  and  that  the  amount  of  upwelled 
radiance  from  any  depth  deeper  than  the  first  attenuation  length  Is  negligible. 

Figure  2  Illustrates  the  Inverse  relationship  of  attenuation  length  with  the 
diffuse  attenuation  coefficient,  Mk.M  Normally,  for  clear  ocean  waters,  the  attenu¬ 
ation  length  Is  20  m  or  greater,  whereas  In  coastal  areas,  the  attenuation  length  Is 
typically  2-4  m.  Although  the  value  of  "k“  depends  on  the  vertical  distribution  of 
scattering  and  absorption  properties,  the  value  of  the  upwelllng  radiation  does  not 
reveal  the  depth  at  which  the  optical  variability  occurs.  In  other  words.  It  Is 
possible  to  have  the  sane  water-leaving  radiance  at  the  surface  for  two  separate 
areas  whose  optical  properties  are  vertically  different  within  the  first  attenuation 
length  (Gordon,  et  al.  1980b).  A  relatively  turbid  homogeneous  water  mass  can  dis¬ 
play  similar  surface  radiance  as  a  clear  water  mass  with  a  turbid  layer  within  the 
attenuation  depth.  Thus,  any  vertical  variation  In  the  optic  properties  that  occurs 
In  the  first  attenuation  length  can  have  a  significant  Influence  on  the  upwelllng 
water-leaving  radiation  sensed  by  the  ocean  color  scanner. 

The  phytoplankton  pigment  concentration  of  ocean  waters  Is  largely  responsible 
for  the  optical  properties.  Ibis  pigment  Is  largely  representative  of  chlorophyll  a 
and  also  Includes  the  phaeoplgment  concentration  (Gordon,  et  al.  1983).  The  spectral 
absorption  characteristics  of  the  pigment  has  significant  effects  on  the  upwelllng 
radiance  at  the  sea  surface  and  consequently  on  the  signal  sensed  by  CZCS.  The 
strong  correlation  between  this  biological  pigment  and  the  water  optical  properties 
have  been  referred  to  as  bio-optical  properties. 

II.  COASTAL  ZONE  COLOR  SCANNER  IMAGERY 

The  Nimbus-7  CZCS  Is  a  6-channel,  multlspectral  scanning  sensor  having  five 
narrow  (+10  nm)  visible  channels  centered  at  443,  520,  550,  670,  and  750  nm  (Hovis, 
et  al.  1978).  The  sensitivity  of  the  spectral  sensors  Is  sufficient  to  measure  sub¬ 
tle  water  color  changes.  The  ground  resolution  at  nadir  Is  800  m,  and  the  satellite 
has  a  sunsynchronous  noontime  orbit  so  that  for  the  Alboran  Sea,  three  consecutive 
days  of  data  collection  are  followed  by  2  days  when  no  data  are  collected.  The 
satellite  acquires  data  at  approximately  the  same  time  each  day,  11:30  hours  local 
solar  time. 

The  quantitative  oceanographic  analyses  of  CZCS  data  require  that  the  data  be 
calibrated  In  units  of  absolute  radiance.  Furthermore,  a  reliable  method  of  separat¬ 
ing  atmospheric  path  radiance  from  ocean  radiance  Is  required.  Gordon  (1978)  and 


Gordon  and  Clark  (1980a)  developed  a  method  that  assumes  the  CZCS  670  nm  channel  can 
be  used  as  a  measure  of  the  atmospheric  aerosol  concentration.  Through  a  weighted 
subtraction  of  this  channel  with  appropriate  subtraction  of  Rayleigh  scattering,  the 
effective  upwelllng  radiance  from  the  ocean  surface  can  be  computed  for  the  CZCS 
443,  520,  and  550  nm  channels.  However,  there  are  several  problems  associated  with 
this  correction  method  that  must  be  understood  In  relation  to  this  study. 

The  first  of  these  problems  Involves  the  proper  selection  of  the  Angstrom 
exponent.  The  Angstrom  exponent  relates  the  wavelength  dependence  of  the  aerosol 
optical  thickness  and  single  scattering  albedo  of  the  aerosols.  (For  a  more  complete 
explanation  see  Gordon  and  Clark,  1980a).  Gordon  and  Clark  (1981)  suggests  that  In 
any  CZCS  Image,  the  radiance  values  of  clear  water  be  used  to  obtain  this  value. 
However,  this  technique  requires  knowing  the  location  of  clear  water  areas  within 
the  CZCS  Image.  Furthermore,  the  assumption  Is  made  that  the  Angstrom  coefficient 
over  the  clear  water  Is  the  same  throughout  the  Image.  If  this  assumption  Is  not 
valid,  the  selected  Angstrom  exponent  value  for  one  part  of  the  Image  will  produce 
erroneous  values  when  applied  to  other  portions  of  the  Image  (Arnone,  1983).  To 
overcome  this  difficulty,  an  Improved  method  has  been  developed  (Arnone  and  Holyer, 
1984)  In  which  selection  of  the  Angstrom  exponent  Is  an  Interactive  analyses  pro¬ 
cedure  that  defines  the  optimum  coefficient  for  the  entire  Image  scene.  This  method 
was  used  for  the  sequential  Images  presented  In  this  study.  It  Is  Interesting  to 
note  that  the  method  showed  that  over  the  relatively  small  area  of  the  Alboran  Sea 
(300  by  200  km),  the  coefficient  does  not  substantially  change. 

A  second  problem  with  this  correction  procedure  Is  the  assumption  that  water 
has  zero  upwelllng  radiance  at  670  nm.  This  Is  a  logical  assumption  for  open  ocean 
and  shelf- type  water  masses  where  minimal  concentration  of  suspended  sediments 
occurs.  However,  significant  upwelled  radiance  at  670  nm  does  occur  In  turbid 
coastal  waters  (Smith  and  Wilson,  1981;  Arnone,  1983a;  and  Arnone,  1983b).  For  these 
coastal  areas,  an  underestimation  of  the  absolute  upwelllng  radiance  In  the  correct 
CZCS  channels  will  result  If  this  method  Is  used.  Fortunately,  the  waters  of  the 
Alboran  Sea  are  relatively  clear  and  only  In  data  from  very  nearshore  waters  (1-3 
km)  should  errors  arise. 

The  final  (and  potentially  most  significant)  problem  with  extracting  quanti¬ 
tative  results  from  CZCS  data  comes  from  uncertainties  In  the  radiance  calibration 
coefficients  of  each  of  the  spectral  channels.  The  sensor  sensitivity  of  the  CZCS 
has  been  decaying  since  launch  with  the  result  that  the  original  calibration 
coefficients  are  no  longer  applicable  (Austin,  1982;  Sun,  1983;  and  Gordon,  et  al., 
1983).  For  the  Alboran  Sea  analyses  presented  here,  an  exponential  decay  of  the 
sensitivity  of  the  CZCS  spectral  channels  has  been  Implemented  that  extends  from  the 
time  of  launch  to  the  time  of  the  Alboran  Sea  data  set  (Goddard  Space  Flight  Center, 
Systems  and  Applied  Science  Corporation,  personal  cotmnunl cation) .  In  addition,  for 
the  short  period  of  the  CZCS  Imagery  presented  here,  the  decay  rate  Is  Insignificant 
when  comparing  Individual  Images. 

In  this  study,  the  ratio  of  the  absolute  upwelllng  radiance  for  the  443-  and 
550-nm  channels  (l.e.,  443/550)  Is  used  to  derive  both  the  diffuse  attenuation 
coefficient  "k"  at  490  nm  (Austin  and  Petzold,  1980),  and  the  phytoplankton  con¬ 
centration  (Gordon  and  Clark,  1980a).  This  method  Implies  an  Inherent  relationship 
between  NkN  and  phytoplankton  concentration.  While  this  relationship  would  be 
expected  In  the  open  ocean  and  shelf  waters,  there  are  limitations  In  coastal  areas 
where  "k"  may  be  heavily  Influenced  by  suspended  sediments  (Arnone,  1983a;  and 
Arnone,  1983b).  As  mentioned  earlier,  the  Alboran  Sea  Is  comparatively  clear 
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although  problems  with  high  sediment  concentrations  could  occur  only  within  1  to  3 
km  of  the  coasts. 


The  CZCS  Imagery  presented  throughout  this  study  Is  displayed  as  the  absolute 
value  of  the  ratio  of  the  443-  to  550-nm  channels.  In  addition,  the  Imagery  has  been 
geometrically  registered  to  a  Mercator  projection  with  a  registration  accuracy  of  +1 
km.  Thus,  reliable  measurements  may  be  made  of  the  temporal  variability  between 
sequential  Images.  As  an  aid  to  both  clarity  and  location,  a  land  mask  of  the  area 
with  the  appropriate  latitude/ longitude  grid  Is  superimposed  on  the  Imagery.  The 
Image  processing  was  done  on  an  Interatlve  digital  satellite  processing  system 
located  at  the  Naval  Ocean  Research  and  Development  Activity. 

III.  SURFACE  DATA  COLLECTION 


Shipboard,  measurements  of  optical  and  biological  parameters  were  collected 
along  a  line  extending  south  from  Marbella,  Spain,  (1 at.  36°28'N,  long.  4°47'W)  to 
the  approximate  center  of  the  Alboran  Gyre  (lat.  35°52'N,  long.  4°35'W)  (Figure  1). 
This  line,  called  the  Marbella  Line,  was  designed  to  transect  the  Coastal,  Atlantic 
Inflow,  and  Alboran  Sea  Gyre  waters. 

The  dates  of  the  various  In  situ  parameters  along  the  transect  during  the 
experiment  are  tabulated  In  Table  1. 


TABLE  1.  Data  Summary  During  Experiment 


October  1982 

Parameters 

4 

5 

6  7 

8 

9  10 

11 

12 

13 

14 

Beam  transmittance 

X 

X 

X 

X 

Conductl vl ty/ temperature 

X 

X 

X 

X 

Spectral  diffuse  attenuation 

coefficient 

X 

X 

X 

XBT 

X 

X 

X 

X 

X 

X 

X 

Chlorophyll -A 

X 

X 

X 

X 

X 

X 

X 

Suspended  sediments 

X 

X 

X 

X 

X 

X 

X 

PRT-5-radlatlon  temperature 

X 

X 

X 

X 

X 

X 

X 

X 

X 

CZCS  coverage 

X  X 

X 

X 

X 

X 

Beam  transmittance  was  measured  with  a  modified  Martek  Instrument  Integrated 
with  a  desk  top  calculator.  The  spectral  response  of  the  Instrument  was  centered  at 
490  nm  (+10  nm) .  The  Instrument  Is  considered  reliable  to  1  percent  transmission.  The 
Instrument  was  lowered  to  a  depth  of  50  m  while  continuously  measuring  percent- 
transmission  and  temperature. 

The  spectral  diffuse  attenuation  coefficient  was  measured  by  an  optical  Instru¬ 
ment  specially  designed  by  Research  Support  Instruments.  The  Instrument  consisted  of 
a  deck  unit  and  a  submersible  unit  both  equipped  with  a  photon  counter  and  a  mono¬ 
chromometer  for  scanning  the  visible  spectrum  between  400-800  nm.  The  submersible 
unit  was  equipped  with  an  upward  and  downward  looking  window,  both  viewed  through 
cosine  collectors  to  measure  upwelllng  and  downwelllng  Irradlance.  The  unit  was 


controlled  by  a  HP -85  calculator  and  operated  In  a  scanning  and  time  averaging  mode. 
Upward  and  downward  spectral  "k"  measurements  were  Integrated  over  the  depth  of  one 
attenuation  length. 

Chlorophyll  and  suspended  sediment  concentration  measurements  were  made  along 
the  Marbella  Line  at  discrete  bio-optical  stations.  Techniques  of  Instrumentation  and 
analyses  were  outlined  In  Garcia,  et  al .  (1983).  Chlorophyll  analyses  Included  both 
continuous  flow  through  fluorometrlc  techniques  (In-vivo)  and  extraction  techniques 
from  In  situ  water  samples  at  discrete  depths.  Suspended  sediment  analyses  were  done 
by  filtration  and  grammetrlc  procedures  on  water  samples  collected  from  depth  and 
continuous  surface  samples. 

The  subsurface  thermal  structure  derived  from  data  analyses  obtained  from  XBTs 
and  thermal  sensors  on  the  optical  Instrumentation  provided  the  basis  for  locating 
the  bio-optical  stations  within  water  masses  and  to  couple  the  surface  measurements 
with  the  subsurface  oceanography. 

IV.  RESULTS  AND  DISCUSSIONS 

A.  Coastal  Zone  Color  Scanner  Imagery 

Atmospherically-corrected  CZCS  Imagery  registered  to  a  Mercator  projection  for 
six  sequential  days  (6,  7,  8,  11,  12,  13  October  1982)  are  presented  In  Figure  3. 

The  Imagery  represents  the  ratio  of  the  443-  to  550-nm  channels  with  the  color  table 
applied  in  the  figure.  From  ratio  values,  the  diffuse  attenuation  coefficients,  and 
phytoplankton  pigment  concentrations  (referred  to  as  Chi  on  the  figure),  are  Illus¬ 
trated  In  this  figure.  (Mote  that  the  color  scales  are  not  a  linear  relationship  of 
"k"  and  pigment  concentrations.)  The  color  scale  represents  higher  turbidity  coastal 
water  masses  as  brown,  shelf-  type  waters  as  yellow  to  green,  and  clear  oceanic 
waters  as  blue. 

In  the  CZCS  Imagery,  the  Alboran  Sea  Gyre  Is  observed  to  have  low  “k"  values 
and  low  phytoplankton  concentration  (blue).  The  gyre  Is  encircled  by  waters  of  ele¬ 
vated  turbidity  and  phytoplankton  concentration  especially  along  Its  northern  and 
eastern  frontal  boundaries.  The  result  Is  a  ribbon  of  significant  color  gradient 
that  is  observed  to  move  eastward  In  the  6  to  8  October  imagery.  (The  changing  color 
Intensity  of  the  front  will  be  discussed  In  the  following  section.)  On  11  October, 
the  frontal  characteristics  appear  significantly  different,  especially  on  the  east¬ 
ern  b'  ^ary  where  the  ribbon  appears  to  have  become  more  diffuse.  The  12  and  13 
October  Images  show  complex  sinuous  features  developing  within  the  frontal  region. 

The  13  October  image  also  shows  the  farthest  eastward  frontal  position,  located 
slightly  to  the  west  of  Cap  Tres  Forcas. 

The  general  water  mass  features  seen  In  the  CZCS  imagery  are  also  observable  in 
the  same  locations  In  the  NOAA-7  Infrared  imagery  for  the  same  period  (Figure  4). 
Analyses  of  the  Infrared  Imagery  during  the  ?Donde  Va?  experiment  by  La  Violette 
(1984),  Indicate  that  a  series  of  submesoscale  cold  water  tongues  were  advected 
clockwise  around  the  gyre.  This  phenomenon  Is  also  observable  in  the  CZCS  Imagery  In 
Figure  3.  Originating  at  the  Strait  of  Gibraltar,  turbid  water  masses  are  shown  to 
enter  coincident  with  cold  water  tongues  that  develop  as  part  of  the  Atlantic 
Inflow.  In  the  Infrared  Imagery  for  the  following  days,  the  tongue  is  advected 
around  the  periphery  of  the  gyre.  In  the  12  October  Image,  the  start  of  one  of  the 
features  Is  readily  observed  as  the  yellow  coloration  apparent  within  the  strait. 

Upon  entering  the  Alboran  Sea,  these  turbid  water  tongues  of  Atlantic  Inflow  water 


36°30'N 


Atmospherically-corrected  and  spectrally  ratioed  (443:550  nm)  sequence 
of  Nimbus-7  CZCS  Imagery  showing  quantitative  biooptical  properties. 
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surface  temperatures  from  NOAA-7  AVHRR 
CS  imagery  shown  in  Figure  3. 


are  not  readily  distinguishable  In  color  from  the  color  of  the  coastal  water,  which 
also  appears  turbid  with  elevated  phytoplankton  concentrations.  As  each  particular 
feature  continues  Its  eastward  migration.  It  maintains  its  Identity  and  sharp 
frontal  boundaries.  Note  that  several  of  these  turbid,  submesoscale  tongues  can  be 
observed  In  a  single  Image.  Their  Initial  speed  appears  to  be  slow  near  Gibraltar 
and  to  Increase  to  an  average  speed  of  approximately  40  cm/ sec  as  they  move  around 
the  gyre.  The  cyclic  appearance  of  these  tongues  appears  to  be  about  12  hours  and  Is 
believed  to  be  related  to  tidal  movement  through  the  Strait  of  Gibraltar  (La 
Vlolette,  1984). 

The  coastal  waters  bio-optical  properties  surrounding  Cap  Tres  Forcas  rapidly 
change  In  the  Figure  3  Imagery.  The  7  and  8  October  Imagery  show  small  westward  flow 
of  the  coastal  waters  as  a  turbid  plume.  On  11  October,  however,  the  plume  developed 
Into  a  thin  ribbon  of  higher  turbidity  water  propagating  northwestward  approximately 
100  km  offshore. 

It  Is  Interesting  to  note  that  In  the  Atlantic,  the  coastal  waters  along  the 
western  Spanish  coast  are  characterl zed  by  high  turbidity  and  phytoplankton  concen¬ 
tration  In  all  six  Images.  The  Rio  Guadalquivir  Is  shown  to  discharge  high  suspended 
sediments  and  nutrients  Into  the  Sea  of  Cadiz.  In  the  sequence  of  Imagery,  this 
turbid  water  mass  appears  to  enter  the  Alboran  along  the  northern  side  of  the  Strait 
of  Gibraltar  and  may  be  a  major  source  of  nutrients  for  the  biological  processes 
taking  place  within  the  Atlantic  Inflow  water. 

The  most  striking  feature  of  the  CZCS  sequence  Is  the  extreme  color  changes 
that  occur  from  1  dayrs  image  to  the  next.  Note  that  color  gradients  define  the 
frontal  positions  of  the  water  masses  In  all  of  the  Imagery.  However,  the  actual 
color  values  In  each  Image  show  strong  changes  that  must  represent  changes  In 
upwelllng  radiation  (assuming  no  errors  Introduced  by  atmospheric  correction).  This, 
In  turn.  Implies  that  optical  or  biological  changes  have  occurred  within  the  first 
attenuation  length.  The  reddish  ribbon  of  Atlantic  Inflow  water  in  the  6  October 
image  changes  to  a  light  blue  ribbon  In  the  7  October  Image,  and  Is  barely  visible 
in  the  8  October  image.  This  indicates  a  day-to-day  change  of  the  diffuse  atten¬ 
uation  coefficient  from  0-15  to  0.07  to  0.06  or  a  phytoplankton  decrease  of  from 
1.05  to  0.45  to  0.35  mg/m  .  A  similar  progression,  though  not  as  dramatic,  may  be 
observed  to  occur  in  the  11  to  12  October  Imagery,  where  “k"  values  decreased  from 
0.09  to  0.075  (little  change  is  found  In  the  13  October  image.).  Notice  that  the 
Atlantic  coast  shows  high  turbidity  In  all  of  the  images. 

A  reevaluation  of  the  atmospheric  correction  process  Indicates  that  an  Inherent 
error  might  be  responsible  for  the  large  temporal  variability.  The  Angstrom  exponent 
selected  for  each  spectral  channel  In  each  Image  could  create  variability  In  the 
absolute  upwelllng  radiance  values.  If  this  were  true,  then  a  relationship  must 
exist  between  the  Angstrom  exponent  and  the  upwelllng  radiance.  To  test  the  validity 
of  this  relationship,  a  pixel  located  approximately  6  km  off  of  Marbella,  Spain,  was 
selected.  The  absolute  upwelllng  radiance  for  the  three  channels  (443,  520,  and  550 
nm)  and  the  associated  Angstrom  exponent  used  for  the  atmospheric  correction  process 
for  each  of  the  six  CZCS  Images  was  plotted  (Figure  5).  The  obvious  lack  of  correla¬ 
tion  In  the  plotted  data  Indicates  the  observed  color  changes  are  not  the  result  of 
the  selection  of  the  Angstrom  exponent.  Thus,  the  Angstrom  exponent  selection  In 
atmospheric  correction  process  Is  not  the  cause  of  the  color  variability  apparent  In 
the  Imagery. 


Notice  that  In  Figure  5  the  upwelllng  radiance  normally  associated  with 
suspended  sediment  concentration  (5S0  nm)  Is  clustered  around  0.038  W/cwr  +.004.  The 
compact  nature  of  this  cluster  suggests  that  the  suspended  sediment  concentration 
had  minimal  changes  at  this  location  In  all  six  Images.  Conversely,  the  large 
changes  In  the  upwelllng  radiance  In  the  443-nm  channel  (0.074-0.032  W/csr)  Indicate 
that  upwelllng  radiance  values  In  this  spectral  Image  were  primarily  responsible  for 
the  color  variability.  Since  the  443-nm  channel  Is  nomally  associated  with  chloro¬ 
phyll  absorption,  the  dally  changes  observed  appear  to  result  from  actual 
chlorophyll  variability. 

The  large  dally  changes  In  the  absolute  upwelllng  spectral  radiance  could  be  a 
result  of  several  factors.  They  could  be  a  result  of  either  the  rapid  bloom  (or 
decay)  of  large  phytoplankton  concentrations  In  the  upper  layers  of  the  ocean,  or 
the  upward  migration  or  upwelllng  advectlon  (or  sinking)  of  the  phytoplankton  popu¬ 
lations.  The  growth  rate  of  phytoplankton  Is  dependent  upon  nutrient  availability, 
temperature,  zooplankton  grazing,  and  Irradlance  distribution.  Cullen  and  Eppley 
(1981)  found  that  the  phytoplankton  concentrations  can  double  per  day  for  a  specific 
vertical  depth.  As  a  result,  phytoplankton  concentrations  can  change  both  horizon¬ 
tally  and  vertically  with  time.  However,  only  those  changes  that  occur  within  the 
first  attenuation  length  affect  the  spectral  upwelllng  radiance  sensed  by  visible 
remote  sensors. 

B.  Ship  Data  along  the  Marbella  Line 

Cross  sections  along  the  Marbella  Line  of  temperature,  beam  transmittance, 
chlorophyll,  and  suspended  sediment  concentration  are  presented  In  Figures  6  and  7 
for  11  and  12  October.  The  temperature  section  on  11  October  (Figure  6)  shows  that  a 
sharp  surface  front  that  sloped  to  the  south  was  present  15  km  from  the  coast.  The 
beam  transmittance  cross  section  for  11  October  closely  resembles  the  temperature 
cross  section.  The  surface  transmittance  front  corresponds  to  the  thermal  front. 
Increasing  from  701  at  the  coast  to  80%  at  the  front.  Further  offshore,  the  water 
became  even  clearer,  with  transmittance  values  of  90+%.  The  vertical  distribution  of 
transmittance  values  shows  an  upward  Intrusion  of  clear  90+%  water  occurred  at  25-m 
depth  25  km  from  shore.  This  clear  water  approached  the  coast  at  approximately  50  m, 
underlying  the  more  turbid  coastal  waters.  Again  the  Isotherms  also  showed  a  similar 
trend  Indicating  strong  Interrelationship. 

The  11  October  chlorophyll  cross  section  In  Figure  6  shows  distributions 
similar  to  that  of  beam  transmittance  and  temperature.  Surface  chlorophyll  concen¬ 
trations  change  at  the  surface  front  from  0.3  to  0.2  mg/1.,  with  higher  chlorophyll 
values  observable  at  depth.  The  strong  chlorophyll  front  corresponds  to  the  90% 
Isollne  of  beam  transmittance  In  the  area  of  the  upward  Intrusion  of  the  clearer 
waters.  In  the  coastal  waters,  the  upper  35  m  contain  high  chlorophyll  concentra¬ 
tions,  while  at  depth  the  exceptionally  clear  water  are  shown  associated  with 
extremely  low  temperature  water. 

The  12  October  temperature  cross  sections  In  Figure  7  show  significant  varia¬ 
tions  from  the  previous  day.  The  surface  thermal  front  has  dissipated  and  a  gradual 
temperature  change  from  17  to  19°C  Is  observed  out  to  50  km  offshore.  The  warm 
(>20°C)  surface  waters  of  the  gyre  are  not  seen  In  this  section,  although  a  gentle 
southern  dipping  of  the  Isotherms  are  observed.  The  change  In  beam  transmittance  for 
the  12  October  cross  section  Is  similar  to  that  of  temperature.  Low  transmittance 
(80%)  measurements  are  observable  40  km  from  the  coast.  The  Isollnes  of  beam  trans¬ 
mittance  like  those  of  temperature  slope  down  to  the  south. 
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Water-leaving  radiance  calculated  from  CZCS  at  a  single  pixel  location 
and  the  Angstrom  coefficient  used  In  atmospheric  correction  for  the 
spectral  channels. 


The  cross  section  of  the  chlorophyll  concentration  for  12  October  also  shows  a 
significant  change  from  the  previous  day.  The  section  shows  exceedingly  high  values 
In  excess  of  2  mg/nr  close  to  the  coast  In  contrast  to  the  0.5  mg/in'3  the3prev1ous 
day.  Surface  waters  40  km  from  the  coast  have  concentrations  of  0.3  mg/nr  compared 
to  0.1  mg/nr  the  previous  d«y.  At  25-m  depth,  a3major  vertical  chlorophyll  concen¬ 
tration  change  occurs.  Maximum  values  of  1  mg/nr  are  Illustrated  to  extend  offshore 
to  40  km  on  12  October.  The  previous  day  this  value  of  1  mg/nr  extends  only  18  km 
offshore.  Notice  that  at  35  m  low  chlorophyll  values  are  observed  on  11  and  12 
October  Indicating  that  temporal  variability  appears  confined  to  the  surface  waters. 

The  changes  In  beam  transmittance,  temperature,  and  chlorophyll  concentration 
for  these  2  days  Indicate  that  an  apparent  southern  movement  of  the  upper  30  m  of 
water  takes  place  In  24  hr.  The  cause  of  this  movement,  as  well  as  the  direction  of 
actual  displacement,  will  be  dlscused  In  the  following  section  when  the  vertical 
sections  of  the  ship  data  are  related  to  the  horizontal  displays  of  the  satellite 
data. 

It  Is  Important  to  note  that  the  vertical  profiles  of  beam  transmittance  pro¬ 
vide  a  clue  to  the  relationship  of  optical  property  variations  In  the  first  attenua¬ 
tion  length  and  the  visible,  remote  sensing  signals. 

Vertical  profiles  of  beam  transmittance  along  the  Marbella  Line  are  Illustrated 
In  Figures  8  and  9  for  11  and  12  October.  The  four  profiles  In  Figure  8  are  repre¬ 
sentative  of  different  water  masses  beginning  near  the  coast  (1)  and  ending  offshore 
In  Alboran  Qyre  water  (4).  The  Alboran  Sea  Qyre  water  has  high  transparency  (97%) 
and  Is  relatively  homogeneous  to  Its  first  attenuation  length  of  21  m.  Profile  (3) 
representing  the  southern  side  of  the  Atlantic  Inflow  water,  has  reduced  transmis¬ 
sivity  (90%)  and  a  nonhomogeneous  upper  layer.  Decreasing  transmissivity  is  observed 
from  the  surface  with  a  slightly  turbid  layer  (88%)  existing  at  12  m.  The  first 
attenuation  length  extended  Into  the  upper  2  m  of  this  layer  such  that  the  water¬ 
leaving  radiation  has  a  spectral  component  that  arose  from  within  this  layer. 

Profile  (2),  located  north  of  the  Atlantic  Inflow,  has  reduced  transmissivity 
(83%)  and  shows  a  significant  variation  within  the  first  attenuation  length.  The 
upper  10  m  Is  optically  homogeneous.  A  sharp  decrease  In  transmissivity  from  85  to 
70%  occurs  between  10-12  m.  The  70%  transmissivity  layer  extends  to  29  m  where  an 
abrupt  Increase  to  92%  transmissivity  occurs  with  depth.  The  first  attenuation 
length  extends  to  17  m  and  Is  partially  Inclusive  of  the  70%  transmissivity  layer. 

The  final  profile  (1)  located  just  offshore  of  Marbella  at  the  start  of  the 
Marbella  Line  has  reduced  transmissivity  In  the  upper  5  m  (67%).  The  Increase  and 
rapid  decrease  In  transmissivity  with  depth  Indicates  a  complex  optical  nature 
within  the  upper  35  m.  However,  since  the  first  attenuation  length  Is  confined  to 
the  upper  8  m,  most  of  this  complex  vertical  structure  does  not  substantially 
Influence  the  water-leaving  radiation  and,  hence.  Is  not  represented  In  the  CZCS 
data. 

Figure  9  representing  profiles  for  12  October  Illustrates  a  profile  located 
within  the  same  coastal  water  mass  (profile  1)  as  the  previous  day.  Within  the  first 
attenuation  length  of  8-m,  surface  transmittance  (65%)  decreases  to  60%.  The  strong 
vertical  variations  have  a  similar  shape  to  the  previous  day.  However,  these  are 
also  outside  the  depth  detectable  by  the  CZCS  sensors.  In  comparing  the  11  and  12 
October  profiles,  note  that  the  vertical  variations  of  the  optical  properties  have  a 
direct  Influence  on  the  depth  to  which  the  CZCS  sensor  can  receive  Its  signal. 
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Figure  9:  Vertical  profiles  along  the  Marbella  Line  of  beam  transmittance  on  12 
October. 


Profile  (2)  In  Figure  9  is  observed  40  km  offshore  In  a  similar  location  as  profile 
(3)  for  11  October  (Figure  8).  Profile  (2)  exhibits  an  optically  homogeneous  Mater 
mass  to  14  m,  which  Is  also  the  first  attenuation  length.  A  lower  transmissivity 
layer  (60%)  Is  observed  between  18  to  32  m.  Through  comparison  with  the  previous 
day,  the  shape  of  profile  (2)  Is  similar  to  waters  found  closer  to  shore. 

From  this  it  can  be  seen  that  the  vertical  variations  of  the  optical  properties 
have  significant  Influence  on  the  upwelllng  radiation  detectable  by  the  CZCS  visible 
range  spectral  sensors.  Since  the  signal  Is  basically  confined  to  the  first  attenua¬ 
tion  length  of  the  water  column,  the  movement  of  the  vertical  optical  layers  Into 
and  out  of  this  attenuation  length  will  result  In  different  measurements  of  upwelled 
radiance.  The  difficulty  Is  further  compounded  since  these  movements  may  change  the 
depth  of  the  first  attenuation  length.  All  of  this  contributes  to  the  high  temporal 
variability  observed  In  the  sequence  of  CZCS  Imagery  In  Figure  3.  An  Important  fact 
to  be  reemphasized  Is  that  changes  below  the  first  attenuation  length  will  not  be 
detected  by  the  visible  range  sensors. 

C.  Comparison  of  Satellite  and  Ship  Data 

The  movement  of  the  coastal  water  noted  In  the  ship  data  is  also  observable  In 
the  11  and  12  October  satellite  data.  The  two  data  sets  show  that  the  complex  verti¬ 
cal  bio-optical  properties  In  the  upper  water  mass  are  likely  to  have  been  responsi¬ 
ble  for  the  high  temporal  and  spatial  variability  observed  In  the  sequential  CZCS 
Imagery.  Most  Important,  the  two  data  sets  add  critical  dimensional  Information  to 
one  another.  For  example,  the  cross  sections  provided  by  the  ship  data  along  the 
Marbella  Line  are  much  more  Informative  when  placed  In  position  In  the  horizontal 
frame  of  the  satellite  data.  Now,  the  apparent  rapid  southerly  movement  of  the 
coastal  water  mass  previously  observed  In  the  ship  data  can  be  seen  to  be  actually 
attributable  to  an  eastward  advectlon  of  the  highly  turbid,  cold  water  feature 
observable  crossing  the  Marbella  Line  In  the  11  and  12  October  satellite  Imagery. 

The  ship  data.  In  turn,  adds  vertical  dimension  to  the  satellite  display  by  showing 
that  the  distribution  of  bio-optical  properties  is  In  the  upper  30  m. 

The  spectral  upwelllng  and  downwelllng  diffuse  attenuation  coefficients  for  11 
and  12  October  are  Illustrated  In  Figures  10  and  11.  These  spectral  measurements 
represent  the  water  mass  variations  across  the  Marbella  Line  coastal  (I)  and  40  km 
offshore  (II)  on  12  October,  Inflowing  Atlantic  water  (I),  and  Alboran  tyre  (II)  on 
II  October.  (The  positions  of  these  stations  are  shown  In  Figures  6  and  7.)  A  shift 
In  wavelength  for  maximum  light  transmission  occurs  In  progressing  north  from  the 
waters  of  the  gyre  to  the  Spanish  coast  (490  to  550  nm)  In  addition  to  an  Increase 
In  the  diffuse  attenuation  coefficient  values. 

The  diffuse  attenuation  coefficient  at  490  nm  (computed  from  the  CZCS  data)  are 
compared  to  ship  observations  In  Figure  12.  There  are  difficulties  In  this  compar¬ 
ison.  First,  the  optical  properties  of  the  regional  water  masses  exhibited  high 
spatial  variability  within  the  800-m  resolution  of  the  CZCS  Imagery.  This  Is  further 
compounded  by  the  1-2  km  registration  accuracy  of  the  satellite  processing  (the  ship 
navigation  also  contained  +2  km  position  errors  when  more  than  10  km  offshore).  In 
addition,  many  envIronmentaT  affects  such  as  sea  state  conditions,  solar  altitude, 
and  ship  shadow  can  contribute  as  much  as  10%  error  to  the  ship  data.  Despite  all  of 
these  factors  the  ship  and  satellite  data  show  good  agreement. 

A  similar  comparison  of  the  phytoplankton  concentration  calculated  from  CZCS 
and  ship  measured  values  of  the  surface  waters  Illustrated  In  Figure  13  also  shows 
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Figure  13:  Comparison  of  ship  measurements  and  Nimbus-7  CZCS  calculated 
phytoplankton  concentrations. 


good  agreement.  Similar  problems  that  Mere  discussed  for  the  MkM  correlation  are 
also  applicable  to  this  figure. 

Profiles  along  the  Marbella  Line  for  the  6-day  sequence  of  the  satellite- 
derived  sea  surface  temperature  and  diffuse  attenuation  coefficient  are  shown  In 
Figure  14.  The  current  velocity  measurements  made  by  the  drift  tracks  of  sonobouys 
(La  Vlolette,  1983)  are  also  Illustrated  for  11  and  13  October.  A  strong  correlation 
between  the  temperature  and  optical  properties  can  be  observed.  Any  small  lack  of 
correlation  between  the  two  satellite  data  sets  can  result  from  several  sources.  The 
most  obvious  of  these  Is  the  1-2  km  registration  error.  Not  as  obvious,  but  equally 
a  source  of  any  displacement  In  the  two  data  sets  Is  the  fact  that  the  two  satellite 
passes  are  2.5  hr  apart.  The  Atlantic  Inflow  water  moves  rapidly  (1.4  m/sec.  Kinder 
1983)  such  that  the  2.5-hr  difference  between  the  CZCS  and  NOAA-7  passes  can  result 
In  a  noticeable  geographic  offset  In  the  values  computed  along  the  line.  Notice  that 
although  the  positions  of  the  surface  thermal  fronts  have  corresponding  optical 
fronts,  the  sea  surface  temperatures  In  the  near  shore  area  does  not  show  the  high 
variability  shown  In  the  optical  properties.  Coastal  and  Atlantic  Inflow  water 
masses  have  similar  surface  temperatures,  yet  some  distribution  can  be  established 
by  the  water  color  signatures.  This  Is  most  apparent  for  the  profile  for  7  October. 

The  positions  of  coastal,  Atlantic  Inflow,  and  Alboran  tyre  water  along  the 
Marbella  Line  are  described  by  the  drift  trajectories  for  11  and  13  October.  The 
current  velocity  observed  on  11  October  Indicates  a  strong  north  component.  These 
current  directions  preceeded  the  leading  edge  of  an  eastward  advectlng  turbid,  cold, 
water  tongue  observable  In  the  satellite  Imagery.  On  this  day  a  strong  eastward 
component  Is  observed  close  to  the  coast  coinciding  with  the  maximum  temperature  and 
optical  front. 

The  13  October  buoy  trajectories  Indicate  that  strong  eastward  velocities  that 
are  located  approximately  30  km  off  the  coast  represent  the  Atlantic  Inflow  (this 
position  coincides  with  the  strong  temperature  and  color  gradients).  It  Is  not 
possible  to  establish  precise  temperature  or  water  color  parameters  that  define  the 
Atlantic  Inflow,  since  the  high  turbulence  Imposed  by  the  variable  Inflow  results  In 
strong  horizontal  and  vertical  mixing  In  the  upper  water  column.  However,  what  does 
appear  consistent  Is  that  the  strong  gradients  observed  In  satellite  Imagery  are 
correlated  with  the  position  of  the  strong  eastward  velocity  measurement. 

V.  SUMMARY  AND  CONCLUSIONS 


The  circulation  of  the  Alboran  Sea  tyre  has  been  characterized  by  a  dally 
sequence  of  visible  (Nimbus-7  CZCS)  and  thermal  IR  (NOAA-7  AVHRR)  satellite  Imagery. 
The  atmospherically-corrected  and  geometrically-registered  Imagery  Illustrates  the 
rapid  north-south  movement  of  the  Atlantic  Inflow  water  as  well  as  associated 
turbid,  cold  water  tongues  that  are  advected  around  the  Alboran  tyre.  The  phyto¬ 
plankton  pigment  concentration  and  the  diffuse  attenuation  coefficient  In  the  upper 
surface  waters  have  been  computed  using  the  ratio  of  two  of  the  CZCS  spectral 
channels.  The  mesoscale  features  characterized  by  cooler  temperatures,  higher 
turbidity,  and  Increased  pigment  concentration  are  shown  to  begin  at  the  Strait  of 
Gibraltar  and  to  move  around  the  tyre.  The  Alboran  tyre  water  Is  characterized  by 
warm  temperatures  and  very  clear  water  with  low  pigment  concentration. 

The  quantitative  bio-optical  properties  measured  from  the  Nimbus-7  CZCS  are  In 
good  agreement  with  simultaneously  collected  ship  data.  The  appreciation  of 
utilizing  quantitative  satellite  data  for  determining  the  spatial  and  temporal 
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variability  of  the  bio-optical  data  Is  demonstrated.  The  limitations  existing  In  the 
data  because  of  physical  contralnts  are  also  shown.  The  6-day  sequence  of  visible 
Imagery  suggests  a  high-daily  variability  of  the  bio-optical  properties  occurring 
along  the  frontal  and  coastal  water  masses.  In  certain  frontal  areas  a  doubling  of 
the  diffuse  attenuation  coefficient  apparently  occurred  within  24  hr. 

Ship  measurements  Indicate  that  the  vertical  bio-optical  distribution  Is 
complex  within  the  upper  30  m,  and  that  strong  vertical  changes  occur  on  a  dally 
basis.  Since  the  satellite  signal  Is  limited  to  the  first  attenuation  length,  strong 
vertical  movements  of  bio-optical  properties  Into  and  out  of  this  first  attenuation 
length  are  reflected  In  the  spectral  satellite  signal,  and  this  vertical  movement 
could  be  the  cause  of  the  changes  In  values  shown  by  the  Imagery. 

Improved  understanding  of  the  variability  of  bio-optical  properties  In  frontal 
areas  Is  required.  The  upwelllng  spectral  radiance  should  be  correlated  with  the 
bio-optical  properties  In  order  to  Improve  our  understanding  of  the  remote  sensing 
of  water  color. 
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